

Center /or 7W6u/ence Research 
Proceedings of the Summer Program 1 994 




2 <27 < 4/ 

N95- 21059 


397 

£ 


Anisotropic structure of homogeneous 
turbulence subjected to uniform rotation 

By C. Cambon 1 , N. N. Mansour 2 AND K- D. Squires 3 


Large-eddy simulation results are used to investigate the development of anisotro- 
pies and the possible transition towards a quasi two-dimensional state in rotating 
turbulence at high Reynolds number. The present study demonstrates the existence 
of two transitions that are identified by two Rossby numbers. The first transition 
marks the onset of anisotropic effects and corresponds to a macro Rossby number 
Ro^ (based on a longitudinal integral length scale) near unity. A second transition 
can be defined in terms of a lower bound of a micro-Rossby number Ro“ also near 
unity (defined in this work as the ratio of the rms fluctuating vorticity to background 
vorticity) and corresponds to a continued development of anisotropy but with an 
increasing emergence of those indicators based on the pure two-dimensional com- 
ponent of the flow, e.g., integral length scales measured along the rotation axis. 
Investigation of the vorticity structure shows that the second transition is also 
characterized by an increasing tendency for alignment between the fluctuating vor- 
ticity vector and the basic angular velocity vector with a preference for corrotative 
vorticity. 


1. Introduction 

Turbulence plays an important role in the performance of many engineering sys- 
tems where the mean flow is dominated by rotation, technologically relevant exam- 
ples include the flow in turbomachines as well as reciprocating engines with swirl 
and tumble. The difficulty in evaluating the performance of these systems is due 
primarily to their complex nature and to the poor performance of turbulence models 
used for prediction of the flow fields encountered in these systems. The principal 
reason for deficiencies in turbulence modeling of rotating flows is in turn related to 
the fact that mean rotation has an indirect and subtle effect on quantities such as 
the turbulence kinetic energy and its dissipation rate, i.e., mean rotation does not 
appear explicitly in the transport equations for either quantity. However, account- 
ing for system rotation is crucial since it alters the nonlinear interactions among 
turbulent scales. Thus, a better understanding of the effects of system rotation on 
turbulent flows is needed to properly guide the development of turbulence models. 

A flow particularly well suited for increasing our fundamental understanding is ro- 
tating homogeneous turbulence since the effect of system rotation can be examined 
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independent of other complicating influences such as mean gradients, body forces 
(other than the Coriolis force), wall effects, or organized structures. A topic central 
to studies of rotating homogeneous turbulence is the generation of anisotropy by 
background rotation and possible transition towards a quasi two-dimensional state. 
Important in this regard is that it can be shown that the initial breaking of isotropy 
arises ultimately from nonlinear interactions which are modified by system rota- 
tion (Cambon Jacquin 1989). Thus, the range of Reynolds and Rossby numbers 
for which anisotropy will develop is not easily predicted. Furthermore, the typical 
measures used to describe anisotropy, e.g., that of the Reynolds stress and length 
scales, require a detailed two-point approach which should also include angular de- 
pendence of the distribution of the energy density in wave space. It is also important 
to emphasize that the Taylor-Proudman theorem cannot predict a transition from a 
three-dimensional to a two-dimensional structure at zero Rossby number, but only 
shows that the 'slow manifold’ (i.e., the manifold which is unaffected by the linear 
‘rapid 7 effects of rotation or inertial waves) is the ‘two-dimensional manifold 7 (see 
Proudman 1916, Cambon, Mansour & Godeferd 1994). 

Insight into the development of anisotropic features in rotating turbulence was 
examined using direct numerical simulation (DNS) and large-eddy simulation (LES) 
by Bardina, Ferziger & Rogallo (1984). They found that anisotropy is primarily re- 
flected by the integral length scales whereas the Reynolds stress tensor remained 
quasi-spherical. Bardina et al did not address, however, the underlying anisotropic 
spectral shape nor the relevant parameter regime for anisotropy development. Sev- 
eral investigators have used both theoretical and experimental approaches (Jacquin, 
Leuchter &: Geoffroy 1989, Cambon & Jacquin 1989, and Jacquin, Leuchter, Cam- 
bon & Mathieu 1990) to investigate anisotropic effects in rotating turbulence, in- 
cluding a non-isotropic spectral description by Cambon (1982) which provided fur- 
ther insight into characterizing the key anisotropy indicators and the relevant pa- 
rameter ranges over which anisotropic effects would develop. These studies showed 
that the most important indicators are expressed using the integral length scales 
with vertical (along the basic rotation axis) separation and vertical or horizon- 
tal velocity components. These previous works also support the conjecture that 
anisotropy in rotating turbulence is triggered when a macro-Rossby number Ro L 
becomes smaller than unity. An optimal collapse of measurements was obtained by 
defining Ro L as 


Ro l 


vgj) 

2f hL\ 3 


( 1 ) 


where £2 3 is the angular velocity of the rotating frame (without loss of generality 
we choose the vertical direction z 3 along the rotation vector and in what follows 
we shall drop the subscript from the definition, i.e., fi = ft 3 ), is the velocity 
fluctuation, and () denotes statistical averaging. The integral length scales are 
defined using the two-point correlation 
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where no summation is implied in the above definition and n is the unit vector in 
the direction referred to by the superscript n. 

Evidence of a second transition in the development of anisotropic features in 
rotating turbulence has come mainly from computational approaches, both DNS 
as well as EDQNM, where high rotation rates are more easily attained than in 
laboratory experiments (e.g., see Mansour, Cambon & Speziale 1991a, b). The sec- 
ond transition corresponds to the lower limit of an intermediate range of Rossby 
numbers in which nonlinear, and therefore the non-isotropic, effects of rotation are 
statistically insignificant. It is conjectured that the second transition occurs when a 
micro- Rossby number Ro x ~ u ( /(2flA) (based on a Taylor microscale A) is smaller 
than 1. For Ro x < 1 non-linear interactions are negligible, e.g., nearly zero deriva- 
tive skewness, and the flow can be accurately predicted using Rapid Distortion 
Theory (RDT). It is important to point out that the RDT solution cannot lead to 
the development of anisotropy, e.g., Speziale et al. (1987) have shown that starting 
with an initially isotropic flow, the flow remains isotropic under very rapid rotation. 
Mansour et al. (1991c) have shown that while rotation will change the anisotropy of 
the Reynolds stress tensor due to the scrambling effects of rotation, it does not lead 
to the generation of new anisotropy. It should also be pointed out that one can pre- 
dict a reduction of the dimensionality, which is part of the anisotropy development 
(Cambon & Jacquin 1989), of rotating turbulence by consideration of resonant tri- 
ads selected at very low Rossby number (Waleffe 1991). However, it is not possible 
to predict the range of Rossby numbers over which the reduction may occur nor the 
measure of the manifold of resonant triads relative to all possible triads using this 
approach. Since the Taylor microscale can be defined from different components in 
anisotropic turbulence, similar to the integral length scales, it is more convenient 
to mark the second transition using a micro- Rossby number defined in terms of the 
vorticity 



where uj, = Note that Ro x and Ro w are proportional since u' ~ u'/X. 

Based upon the summary presented above, there exists evidence to support the 
simple scheme of an intermediate range of Rossby numbers over which anisotropies 
in rotating turbulence will develop. This range is marked by a macro- Rossby number 
close to 1 (where the macro- Rossby number is based on “vertical” quantities) for the 
upper bound and a micro-Rossby close to 1 for the lower bound. The intermediate 
range of Rossby numbers involves the Reynolds number through the ratio L/X 
and is also consistent with the classical use of an ‘Oszmidov scale’ (Jacquin et 
al. 1989). The Oszmidov scale for rotation can be defined by Lo = (e/(2D) 3 ) 1/ ' 2 
(see also Zeman 1994) where the Brunt-Waisala frequency has been replaced by 
the basic vorticity 2fl so that Lo characterizes the smallest scale up to which the 
rotation is acting. A value of the macro- Rossby number nearly 1 corresponds to an 
Oszmidov scale close to the integral length scale Lo = L, whereas the micro-Rossby 
number close to unity corresponds to an Oszmidov scale near the Kolmogorov scale 
Lo = r). Hence, the intermediate range of Rossby numbers can be linked to the 
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largest ( L ~ k 3 ? 2 /e) and the smallest length scales (77 = (i/ 3 /e) 1 ^ 4 ) having physical 
relevance in homogeneous turbulence. 

While the majority of previous work has provided evidence of two transitions in 
the development of anisotropy in rotating turbulence, there exists no previous work 
in which both transitions are evident in the same experimental case or computational 
run. It is difficult in experiments to obtain large enough rotation rates to observe 
the second transition while previous computations have not had statistical samples 
large enough to observe the first transition. An additional difficulty associated with 
DNS is that the Reynolds number range achievable in DNS prohibits a detailed 
examination of the generation of anisotropy arising from non-linear interactions 
since molecular viscosity effectively drives the flow to the RDT limit before non- 
linear effects become significant. Thus, the most viable approach for examining 
anisotropy development is simulations performed at high Reynolds numbers using 
large-eddy simulation. Relevant in this regard is the recent LES study by Bartello, 
Metais and Lesieur (1994). Bartello et al. examined the structural evolution of 
initially isotropic turbulence subjected to uniform system rotation and observed 
a single transition marking the generation of anisotropy. The transition in their 
work occurred around RcP ~ 1 (referred to as the second transition above) with 
significant anisotropy developing only for Ro w < 1. 

Therefore, the primary motivation of the work presented in this paper has been 
to provide a detailed examination of the evolution of statistical and structural 
anisotropy indicators in rotating homogeneous turbulence at high Reynolds num- 
bers. Especially important to this work are simulations of rotating homogeneous 
flows in which the intermediate range of Rossby numbers is achieved, i.e., Ro L < 1 
and Ro“ > 1 in order to permit observation of all possible transitions. Of further 
interest in this work is examination of the underlying structure of the vorticity field 
as anisotropy develops. Evidence exists of a tendency for alignment of the vorticity 
along the rotation axis with a preference for alignment of corrotative (i.e., parallel 
or cyclonic) vorticity compared to contrarotative (i.e., antiparallel or anticyclonic) 
vortices (e.g., see Bartello et ai 1994). 

An overview of the LES database used in this work is presented in §2, including 
time development of measures of the anisotropy and examination of the vorticity 
structure using a joint pdf of the cosine of the angle between the X 3 -component 
of the vorticity vector and background rotation vector with the local micro- Rossby 
number u>iUi/(2Q). Discussed in §3 is the anisotropy ‘shape’, i.e., interpretation of 
the ‘cigar type’ of covariance matrices for both the velocity and the vorticity field, 
using extensively the concepts of ‘directional dependence’ and ‘polarization’ (Cam- 
bon, Jacquin &; Lubrano 1992) or, almost equivalently, those of ‘dimensionality’ 
and ‘componentality’ (Reynolds Sz Kassinos 1994). The relevance of a classic cen- 
trifugal instability for predicting the asymmetry of the vorticity structure in terms 
of corrotative and contrarotative vortices is discussed, together with a summary of 
the work from this study, in §4. 
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2. Database overview and simulation results 

2.1. LES database 

In this work a pseudo-spectral method was used to integrate the filtered Navier- 
Stokes equations for an incompressible fluid in a rotating reference frame. The 
formulation uses an implicit treatment of linear rotation terms which allows exact 
integration of the rapid distortion equations (Rogallo 1981, Mansour et al 1991a). 
Subgrid-scale stresses were parameterized using a spectral-eddy viscosity (Kraich- 
nan 1976, Chollet k Lesieur 1981) modified for rotating turbulence (see Cambon, 
Bertoglio k Jeandel 1981 and Squires, Chasnov, Mansour k Cambon 1994). Use of 
a subgrid-scale model yields a third relevant Rossby number 


Ro" = 


8fl 2 


( 4 ) 


referred to as the cut-off Rossby in this manuscript ( k c is the cut-off wavenumber and 
E(k,t) the classic shell- averaged energy spectrum). Since the molecular viscosity 
was dropped from the Navier-Stokes equations, the simulations can be considered 
to represent asymptotically the dynamics at infinite Reynolds number. 

The computations were performed using an elongated domain in which the verti- 
cal dimension is four times the horizontal length. This geometry and the isotropic 
initial data were carefully optimized in order to permit the development of the 
flow without disturbing statistical symmetries. This is especially important in the 
present work since the initial three-dimensional isotropy is broken only by nonlin- 
ear dynamics. Care was also exercised in order to insure that the length scales 
remain small compared to the size of the computational box. Four datasets were 
used in the course of this work: the first case (referred to as DO) was obtained at 
a resolution of 64 x 64 x 256 collocation points and rotation rate of 0.1. Case DO 
was used primarily as a means of obtaining preliminary structural results but also 
permits examination of the effect of resolution. The other cases, referred to as Dl, 
D2, and D3 for computations performed using rotation rates of 0.1, 0.5, and 1.0, 
respectively, were performed using resolutions of 128 x 128 x 512 collocation points. 
The reader is referred to Squires, Chasnov, Mansour k Cambon (1993, 1994) for 
additional details concerning the simulations. 

2.2. Statistical evolution 

Shown in Fig. 1 are the time histories of the three relevant Rossby numbers Ro L , 
JZo^, and Ro 99 for Cases DO and Dl (Fig. la), Case D2 (Fig. 16), and Case D3 
(Fig lc). The time axis in Fig. 1 has been made dimensionless by the initial eddy 
turnover time. It is evident from the figure, especially for Cases DO and Dl, that 
the initial value of any Rossby number is large enough to ensure that the rotation 
has an initially weak influence and that there is adequate time for development of 
triple correlations as in non- rotating turbulence. We find that els in previous DNS 
and experimental studies, the value of the macro- and the micro-Rossby numbers 
decrease with time. Eventually, the intermediate range of Rossby numbers, i.e., 
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FIGURE 1a. Rossby number development, (a) Case DO: « © Ro L \ Ro w ; 

Ro* 9 . Case Dl: Ro L \ Ro u ; Ro* 9 . 

Ro L < 1 and Ro u > 1, is attained during the simulation. Fig. la shows that the 
intermediate range is maintained for more than a decade for Case Dl, ending at 
approximately the final time of the calculation. In addition, the results in Fig. 1 
also show the effect of resolution (Cases DO and Dl) as well as a portion of the 
calculation occurring after the second transition (Cases D2 and D3). 

To test the validity of computing statistics using only the large scales in an LES 
computation, Chasnov (1994) assumed that an inertial range extends towards infin- 
ity after the cut-off wavenumber. Chasnov found that this extension of the spectrum 
yields relatively small corrections to the values of the kinetic energy and the integral 
length scales as compared to computing these values from the resolved scales. Thus, 
computing the macro- Rossby number using data from only the resolved scales does 
not present any theoretical or practical difficulties. Because the smallest scales of 
motion are not represented in large-eddy simulation, however, the definition of a 
micro-Rossby deserves additional discussion. It is worth noting that an extension 
of the spectrum past the cut-off wavenumber to account for the unresolved scales 
will result in a divergent vorticity field; therefore, this approach cannot be used to 
define a micro-Rossby number. Thus, in this work the micro-Rossby number shall 
be defined using the vorticity field from the resolved scales, yielding an effective Ro w 
in the same way as the use of an eddy viscosity in LES yields an effective Reynolds 
number. 

Note that it is also possible to define a pseudo- Reynolds number Re* based on the 
resolved scales as the ratio between the macro- and the micro-Rossby numbers in 
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agreement with the isotropic relationship, i.e., Re* = k 2 /(ue) = (2/3)(i?o u, / Ro L ) 2 . 
The time development of Re* is presented in Fig. 2. As opposed to the behavior 
observed in DNS and experiments, Re* increases with elapsed time since the integral 
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*/r( 0) 

FIGURE 2. Time development of “Reynolds number”. Case Dl; Case 

D2; Case D3. 


length scale continually increases while the resolved Taylor microscale (or A oc u'/lo') 
tends to be blocked by the cut-off scale 2 n/k c . The latter behavior is consistent 
with the nearly equivalent values of Ro“ and Ro * 9 , with the exception of very early 
times in the development (Fig. 1). Some discrepancies between Cases DO and Dl 
shown in Fig. la illustrate the effect of the resolution: Ro w and Ro S9 axe shown to 
be nearly equal during the entire time development only at the finer resolutions. 

Of the various statistical measures of anisotropy, the measure most widely uti- 
lized in previous investigations has been that of the Reynolds stress, i.e., bij = 
Rij/Ru ~ $ij/ 3, where R l} = Note that there is only one component rel- 

evant in the semi-axisymmetric configuration created by rotation, 633 = bjj ^ 
since 61 1 = 622 = "^33 /2 and the off-diagonal components are zero. Fig. 3 shows 
the evolution of the Reynolds stress anisotropy as a function of time for the three 
high resolution cases considered. We find a mild positive maximum value of about 
0.06 which is consistent with a ‘cigar-type’ shape (i.e., elongation in one direction) 
for the Reynolds stress tensor with a ratio of horizontal and vertical principal com- 
ponents £((uf) + {u\))/{ul) about 0. 7-0.8. A very good collapse of the 3 cases (Dl, 
D2, D3) is demonstrated in Fig. 4 in terms of the macro-Rossby number and the 
anisotropy does indeed appear to be triggered around Ro L = 1. A similar collapse, 
and further confirmation of the first transition, is found using the vorticity ratio 
5 ( (wj ) + ( u - , 2 ) )/ ( a ' , 3 ) ■ ki order to provide an anisotropy measure for the vorticity 
in agreement with the one used for the velocity, the component 633 of the vorticity 
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FIGURE 3. Time development of £>33 component of Reynolds stress anisotropy. 
Case Dl; Case D2; Case D3. 



Figure 4. Development of 633 as a function of the macro Rossby number. 
Case Dl; Case D2; Case D3. 


deviator 


u> 




<^ 3 > 1 
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is plotted in Fig. 5. We find that the vorticity tensor exhibits the same type of 
anisotropy shape (‘cigar- type’) as the Reynolds stress tensor. 



FIGURE 5. Development of vorticity anisotropy 633 as a function of the macro 
Rossby number. Case Dl; Case D2; Case D3. 

Similar anisotropy is also found using the most relevant integral length scales, 
L33 and Xji + L 22 = (Fig. 6). An anisotropy measure of the integral length 
scales can be defined as 


b L 


WHi + K 1 


( 6 ) 


based on a true tensor of ‘2D-energy components’ £fj = (uiUj)L^ (Cambon 1990). 
The quantity b L is shown in Fig. 7 and is also consistent with the results observed 
in Figs. 4-6. 

Finally, it is interesting to note that the slight plateau exhibited in Fig. 3 and 
Fig. 4 is consistent with the blocking of nonlinear anisotropic effects previously 
observed in DNS and EDQNM computations after the second transition Ro u ~ 1. 
However, the fact that the anisotropy reflected by the vorticity ratio and the integral 
length scale ratio continues to increase after the second transition (Ro w < 1) is 
surprising and requires further investigation. 

2.S. Vorticity structure 

The anisotropy reflected by the Reynolds stress tensor and the vorticity corre- 
lation tensor has motivated a refined structural analysis. The joint-probability of 
cos a — u> 3 /u; and of the local Rossby number, y/uiUi/2Sl, is plotted for different 
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FIGURE 6. Development of integral length scales L33 and Ljj as a function 

of the macro Ftossby number. Case Dl: e » L\ x + L\ 2 \ £ 33 * Case D2: 

L\ j + L 22 » ^ 33 * Case D3: L\^ + ^22» -^ 33 * 



FIGURE 7 . Development of length scale anisotropy b L as a function of the macro 
Rossby number. Case Dl; Case D2; Case D3. 
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FIGURE 8. Joint pdf of local Rossby number and cosine of the angle between 
vertical vorticity and vorticity magnitude for Case D2, t/r(0) = 13. Contour levels 
drawn between 0 and 8 with an interval of 0.4. 



FIGURE 9. Joint pdf of local Rossby number and cosine of the angle between 
vertical vorticity and vorticity magnitude for Case D2, t/r( 0) = 48. Contour levels 
drawn between 0 and 8 with an interval of 0.4. 

times in Figs. 8 through 11. The x-axis represents binning in terms of the cosine 
of the angle between the rotation axis and the vorticity vector. Values of the co- 
sine equal to —1 represent contrarotative or antiparallel eddies, values of the cosine 
equal to 1 represent corrotative or parallel eddies. The use of the cosine of the angle 
rather than the angle itself ensures equal probability of occurrence in any bin. The 
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FIGURE 10. Joint pdf of local Rossby number and cosine of the angle between 
vertical vorticity and vorticity magnitude for Case D2, t/r(0) = 109. Contour levels 
drawn between 0 and 8 with an interval of 0.4. 



FIGURE 11. Joint pdf of local Rossby number and cosine of the angle between 
vertical vorticity and vorticity magnitude for Case D2, t/r(0) = 1168. Contour 
levels drawn between 0 and 8 with an interval of 0.4. 

y-axis represents binning in terms of the local micro- Rossby number, and because 
fi is uniform, it is also a measure of the magnitude of the vorticity fluctuation. 
Fig. 8 corresponds to the beginning of the intermediate range ( Ro L = l,f?o w = 8, 
see Fig. la) and the horizontal layering reflects a complete isotropy at a very fine 
level (statistical equidistribution of the orientation of the vortices). At smaller Ref 
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ranging from 2 to 1 (Figs. 9 to 11 ), the distribution shows a distinct tendency for 
the vortices to be increasingly aligned with the basic rotation axis. The distribution 
exhibits two peaks near cos a = 1 and cos a = — 1 with a preference of corrotative 
eddies (cos a = 1 ) clearly exhibited in the distribution. This finding is consistent 
with that previously observed by Bartello et al (1994) that corotating vorticity 
tends to become the dominant vorticity in the flow. 

3. Interpretation of the anisotropic structure 

3.1 . Velocity and vorticity fields 

The development of anisotropy can be easily argued for most homogeneous flows 
based upon the type of mean deformation. For instance, the anisotropy created by a 
distorting duct can be interpreted as follows: in the direction of maximum stretching 
(chosen as X 3 for sake of discussion), the vorticity is elongated so that the intensity 
of vortices is reinforced in accordance with conservation of angular momentum and 
(utj) becomes dominant with respect to other (horizontal) principal components. 
In this case the covariance assumes a ‘cigar type’ shape under the influence 

of the mean strain. Because of the more rapid rotation of the eddy in the plane 
normal to the axis of stretching, the velocity components orthogonal to the X 3 axis 
will become larger, creating a ‘pancake type’ shape for the tensor i.e., two 

dominant components relative to a third weak component. The interesting feature 
of the rotating flow, however, is that both the single-point covariance matrices of 
the velocity and vorticity field possess a ‘cigar type’ shape. This is a very specific 
feature not observed in other anisotropic flows subjected to mean deformation and 
in turn motivates a more refined interpretation than can be obtained from Figs. 3-5. 

The specific anisotropy created by rotation through non-linear interactions can 
best be described using a two-point approach, or at least a single-point approach 
that carries information about both the ‘dimensionality’ and ‘componentality’ of 
the flow. For this purpose it is convenient to utilize spectral variables (e, Z, h ) that 
completely characterize the double correlations at two points for any homogeneous 
anisotropic field. These three variables are invariants of the flow: the energy, e, the 
polarization anisotropy, Z, and the helicity, /i, and correspond to double correlations 
of the two components (or helical modes intensities) in the plane normal to the 
wave vector k (see Cambon & Jacquin 1989 for a complete discussion). In the semi- 
axisymmetric configuration created by rotation through nonlinear effects, only e and 
Z are relevant and they depend only on the wavevector modulus k and on the cosine 
of the angle between k and O, i.e., cos 0 * = k^/k. It is important to note that e is 
a three-dimensional energy spectrum (e cx ~{ 2 * 2 ;)) but is not spherically averaged. 
Thus, the difference between e and the spherically averaged value, i.e., e(&, cos#*) — 
E(k)/(4irk 2 ), characterizes a directional anisotropy. The complex scalar Z gives 
at fixed k the difference of the two nonzero principal components of (u*Uj) as its 
modulus, |Z|, and the orientation of the related principal axes as its argument. 
Thus, Z characterizes a polarization anisotropy. The transport of these quantities 
is governed by 

(- + 2 uk 2 ) e(*,Jfcj/M) = T e (k,k 3 /k,t) 


(7a) 
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(±+ 2 uk 2 + tin^)Z(k,k 3 /k,t) = T z (k,k 3 /k,t) (76) 

which describes the dynamics in the presence of rotation with an exact separation of 
contributions from linear (or rapid) terms on the left hand sides and of contributions 
from nonlinear terms on the right hand sides. The terms T and T are spectral 
transfer terms that involve cubic correlations. In an LES calculation T* and T 
may be replaced by (T e - 2u ag k 2 e) and (T* - 2v ag k 2 Z), so that v is replaced by v, g 
in (7), whereas denote the resolved transfer terms. 

One of the primary advantages of the decomposition into the invariants e and 
\Z\ is that the spectral tensor of double correlations, vorticity as well as velocity, at 
two points can be expressed in terms of t and Z in a fixed frame of reference. Any 
single-point tensor obtained by subsequently integrating in wave-space can then 
be split into a pure (three-dimensional) isotropic contribution (from E/(47rfc 2 )), a 
contribution from the directional anisotropy (e — E/(4nk 2 )), and a contribution 
from the polarization anisotropy Z . For convenience this splitting is applied to the 
Reynolds stress tensor, 






and to the tensors introduced by Reynolds & Kassinos (1994), i.e., the dimension- 
ality structure tensor 


and the circulicity tensor 


D <1 = (t - 2 b h ) 

*w(f- +^-‘«)- 


The superscripts e and 2 in (8)-(10) denote contributions from e - E/{4irk 2 ) and 
from Z, respectively. The dimensionality tensor 


V 2 i/>» = -Wj = ( lla ) 


is equivalently defined by 

D,, = J 2^e(k,t)d'k (1U) 

and displays the directional dependence of the spectral distribution of energy. It is 
clear from (9) and (11) that ‘directional dependence’ and ‘dimensionality’ are similar 
concepts. Eq. (8), however, shows that the ‘componentality’ of the Reynolds stress 
tensor (denoted by its deviatoric part 6 tJ ) results from additive contributions from 
‘dimensionality’ (or ‘directivity’) 6®j and from ‘polarization 6*^. The circulicity 
tensor Fij = (V’i.tV’j.fc) characterizes the contribution from the largest scales to 
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FIGURE 12. Reynolds stress anisotropy for Case Dl, intermediate range of Rossby 
numbers Ro L and i?o w occurring for 100 < t/r( 0) < 3000 (marked by solid vertical 
lines on figure). fe 33 ; 6| 3 ; fe| 3 . 


the vorticity covariance matrix since it is obtained by integrating (Q*Uj)/k 2 (see 
Eq. 11a). It is also worth noting that the ‘directional dependence’ of Fij is exactly 
the same as that of the Reynolds stress tensor whereas its ‘polarization’ is exactly 
the opposite. 
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FIGURE 13. Reynolds stress anisotropy for Case D2, intermediate range of Rossby 
numbers Ro L and Ro w occurring for 15 < </r(0) < 300 (marked by solid vertical 
lines on figure). 633; 633; &33- 


Shown in Fig. 12a (linear time axis) and Fig. 126 (logarithmic time axis) are 
the relevant anisotropy indicators in semi-axisymmetric turbulence, i.e., 633 and 
components 633 and 633, computed using (8), (9), and (116) for Case Dl. The 
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strong dominance of 633 relative to 633 demonstrates that the ‘slow’ effect of rota- 
tion is mainly a ‘dimensional’ one. This is further supported by the observation 
that the spectral density of energy tends to concentrate in the wave-plane k 3 — 0, 
i.e., the pure two-dimensional manifold, since a complete concentration at k 3 = 0 
corresponds to a vanishing of vertical variability in physical space (see Squires et 
al 1993). If a pure two- dimensional/ three- component (2D-3C) state was reached, 
^33 /q 2 would be zero according to (11) so that 6 33 ~ b% 3 would be equal to 1/6. 

It is also important to remark that the contribution from the ‘polarization’ indi- 
cated by a slightly negative value for 6 33 in Fig. 12 is not insignificant. For example, 
the RDT solution, through scrambling effects by inertial waves, modifies the angle 
of polarization through a (fc 3 /k)- dependent angle of rotation (2 (ltk 3 /k) so that the 
contribution from Z is ‘killed off’ on average by linear effects and, hence, b 33 ~ 0. 
However, inertial waves cannot affect Z in the slow manifold k 3 = 0 where Z is real 
and the polarization discriminates contributions to vertical velocity components 
e + Z from contributions to horizontal velocity components e — Z. The quantity 
c — Z at fc 3 = 0 is precisely the two-dimensional/two-component mode since it con- 
tributes only to the horizontal velocity field in the slow (two-dimensional) manifold 
and & 2 (e — Z) contributes only to the vertical vorticity component. Accordingly 
e + Z at fc 3 = 0 is the 2D-3C mode in the slow manifold, which only contributes to 
the vertical velocity component, whereas k 2 (e + Z) contributes only to horizontal 
vorticity components. The dynamics of these two modes of the slow manifold are ac- 
counted for in physical space through products of integral length scales and related 
Reynolds stress components (e.g., Eq. 6) but are weakly reflected by single-point 
velocity or vorticity correlations obtained by averaging over all angles in spectral 
space. In any case, examination of contributions from only the slow manifold would 
indicate that the weak negative value of b 33 may characterize a tendency of the 
velocity field to be rather horizontal and of the vorticity field to be rather vertical. 
For the vorticity correlations the ‘cigar type’ shape created by the ‘dimensionality’ 
for the circulicity tensor is slightly increased by polarization (since b\ 3 - 6| 3 > 6§ 3 ) 
according to (9), whereas the ‘cigar type’ shape of the Reynolds stress tensor is 
slightly diminished since b% 3 + 6| 3 < b e 33 , However, it is important to also note that 
the ‘dimensionality effect’ is not balanced by the ‘polarization effect’. In addition, 
the latter effect can be interpreted as a ‘componentality effect’ but restricted to the 
slow manifold where inertial waves cannot destroy it. 

It is further interesting to examine 6 33 and its contributions from e and Z for 
Cases D2 and D3. The development shown in Figs. 13 and 14 are consistent with 
the previous discussion for Fig. 12 but also clearly suggest that the only significant 
contribution from Z is obtained after the second transition in accordance with the 
location of i?o w = 1 (cf. Fig. 16 and Fig. lc). The stabilization of 6 33 at early times 
following the second transition is due to the balanced increase of 6 33 and decrease 
^33* For increasing time past the second transition it is also clear from Figs. 13 and 
14 that the behavior of 6 33 is driven by the increasingly negative value of b 33 . It 
is possible that for longer evolution times 6 33 may become negative, indicating a 
two- dimensional/ two- component state of the Reynolds stress tensor. However, it is 
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FIGURE 14. Reynolds stress anisotropy for Case D3, intermediate range of Rossby 
numbers Ro L and Ro " occurring for 7 < t/r( 0) < 100 (marked by solid vertical 
lines on figure). 633 ; — ^ 33 » ^ 33 * 


difficult to confirm this tendency even in high resolution LES since the imposition 
of periodic boundary conditions will have an adverse effect on the overall evolution 
of the flow. 

Finally, the three cases are plotted in Fig. 15 in terms of Ro L . The collapse of the 
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data for the three cases is remarkable and demonstrates that the findings discussed 
for Case D1 apply to the other two cases as well. Furthermore, the results in 
Fig. 15 also show that the macro Rossby number is the key parameter determining 
the state of the turbulence. The onset of anisotropy is clearly indicated for Ro L < 1 
while a second transition occurs at approximately Ro L <0.1 which corresponds to 
Ro w <1. 

Analogous quantities to the anisotropy measures of the velocity field shown in 
Fig. 15 have been plotted in Fig. 16 for the vorticity field. ^ These quantities were 

calculated using a k 2 weighting factor so that 63* — 633 yields 633, i.e., the 

anisotropy of the single-point vorticity correlation tensor. The relative role of Z 
(or polarization) in the development of the anisotropy is weak compared to the 
contribution due to the role of e (i.e., dimensionality). This is a reflection of the 
concentration of energy around £3 = 0. As was also observed for the velocity, the 
collapse of the vorticity anisotropy for the three cases in terms of the macro Rossby 
number is again remarkable. In this case the first transition is clear while the second 

° (k 4 * * 7 z) 

is marked by a change in the slope of i> 33 

Finally, it is again worth noting that the ‘cigar-type’ shape created by the pure 

‘dimensionality’ effect ( b e 33 > 0, b[ 3 > 0), has nothing to do with a tendency of 

the velocity vector or the vorticity vector to be aligned with the rotation vector. 
The relative concentration of spectral energy at £3 = 0 corresponds to a lack of 
energy about “vertical” wavenumbers k (i.e., nearly aligned with the rotation axis). 
The contribution of these modes at “vertical” wavenumbers is to horizontal velocity 
or vorticity components since u and Q are normal to k. In other words, 633 > 0 

(and 633 e) > 0) characterize the lessening of dimensionality in the presence of the 
incompressibility constraint. 


4. Summary 

The relevance of two threshold Rossby numbers which delineate anisotropic ef- 

fects in rotating turbulence has been demonstrated using the database from high- 
resolution LES at asymptotically high Reynolds numbers. Anisotropy created by 
background rotation can only be initiated by nonlinear interactions and is triggered 
when a macro- Rossby number becomes smaller than 1. An excellent collapse of the 
main anisotropy indicators was demonstrated in terms of the macro- Rossby number. 

The second transition corresponding to a micro-Rossby number iZo w close to 1 
was found in previous experimental, EDQNM, and DNS results as the beginning of 
the shut-off of the energy cascade indicating that non-isotropic, non-linear effects 
were frozen. In the LES results used in this study, a second transition is indeed 
observed for Ro P < 1 but represents a reorganization of the anisotropy rather than 
a stabilization. The description of this complex anisotropic structure has motivated 
a detailed analysis of the ‘dimensionality’ and ‘polarization’ effects of steady system 
rotation. These effects can be summarized by considering a rectangular domain in 
spectral space where |Ar 3 /Jb| is plotted on the vertical axis and k on the horizon- 
tal axis (Fig. 17). The domain can be separated horizontally in a narrow band 
\h/k\ < e s << 1 near Jt 3 = 0 (the slow manifold denoted ‘s’) and the remaining 



418 


C. Cambon, N. N. Mansour & K. D. Squires 


w 

°L 

D w 

U 

D w 

(0 

.8 

d l 

D s 

u 

D 8 

0) 


k 


FIGURE 17. Sub-domains in fc-space. 


upper band (denoted ‘w’) where the scrambling effect of inertial waves are present. 
The spectral domain can be further subdivided into three bands of wavenumbers 
that characterize the largest scales (L), the energetic scales (u), and the vorticity 
scales (u>). The six sub-domains corresponding to these divisions are denoted as 

in Fig. 17. The first transition is characterized by the 
beginning of a lessening of dimensionality and concentration of energy towards the 
slow manifold in agreement with 633 > 0 (which characterizes an angular drain 
of energy from D™ to D *) and b^ k ^ > 0 (which characterizes an angular drain 
of energy from D™ to D * ). Polarization, or componentality, only affects the slow 
manifold and is reflected by b L > 0 (which characterizes D 3 L ). When the Rossby 
number is in the intermediate range ( Ro l < 1 and Ro * > 1 ), the polarization 

indicators 633 and 633 ; remain weak so that the global anisotropy reflected by 
633 and 633 is an almost pure dimensional effect. After the second transition, the 
polarization effect in the slow manifold seems to emerge not only in D* L , but also in 
D * through 6| 3 < 0 so that 633 = 633 + 633 results from opposite ‘dimensional’ and 
‘polarization’ effects. Nevertheless, the vorticity anisotropy 633 = &! 


— /|(* 2 e) *) 


con- 


tinues to be essentially due to the dimensional contribution 633 e K Such a pattern 
may correspond to opposite values of the polarization anisotropy in the domains 
D% (Z < 0 b L > 0 ), D* u (Z < 0 -> &f 3 > 0) and D* w ( Z > 0). 

An analysis of the vorticity structure was performed by means of a joint-pdf repre- 
sentation of the cos a — and the local Rossby number >/u^u?7/(2Q). The second 
transition (/?o w ~ 1 ) is shown to be characterized by a tendency of the vortices to 
be aligned with the basic rotation vector with an increasing dominance of corrota- 
tive eddies. This last result could be explained by a classical Bradshaw- Richardson 
stability criterion (Bradshaw 1969), also consistent with Rayleigh centrifugal crite- 
rion, at least for values close to zero absolute vorticity u>i + 2 — 0 or zero tilting 
vorticity + 2Qi, as discussed by Cambon et al , (1994a). These results for the 
vorticity structure are also in qualitative agreement with those of Bartello et al. 
(1994). However, in the present study this is a side effect and not the main source 
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of the anisotropy. The fact that the maximum value of u> 3 /(2 fl) is nearly the same 
as Ro w is consistent with the fact that the componentality effect (with corresponds 
to a polarization effect in the slow manifold) has a relatively weak effect on the 
structure of the vorticity field. Statistically, the ‘cigar type’ shape of the vorticity 
correlations tensor mainly reflects the pure dimensional effect. 

It should also be noted that anisotropic trends would be blocked for Ro w < 1, 
except for quantities only dependent on the slow manifold (e.g., 5^), in accordance 
with a shut-off of the energy cascade if viscous effects were present in the simulations. 
The reorganization of the anisotropy (especially as reflected by &| 3 in the Reynolds 
stress as well as in the circulicity tensor) after the second transition is only allowed 
in the absence of dissipation. 
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